Characteristic Energy of the Coulomb Interactions and the Pileup of States 
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Tunneling data on Lai.28Sri.72Mn207 crystals confirm Coulomb interaction effects through 
the dependence of the density of states. Importantly, the data and analysis at high energy, E, 
show a pileup of states: most of the states removed from near the Fermi level are found between 
^ 40 and 130 meV, from which we infer the possibility of universal behavior. The agreement of our 
tunneling data with recent photoemission results further confirms our analysis. 
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Coulomb interactions (CPs) between electrons can 
play a dominant role in strongly correlated systems and 
metal-to-insulation transitions. These include Mott in- 
sulators, complex oxides [H, Q, colossal magnet oresis- 
tive materials y, and other novel metallic materials like 
graphenefsj. Coulomb interactions have also been stud- 
ied as a possible mechanism of the pseudogapjl, 0, Q in 
the density-of-states (DOS) of high-temperature super- 
conductors. The Coulomb repulsion pushes some states 
near Fermi level, /ii?, to higher energies, as described by 
Al'tshuler and Aronov (AA) Despite years of study, 
however, the essential implication of the theory, the con- 
servation of states, has never been shown. In supercon- 
ductivity most of the states depleted from /i^? pile up just 
above the energy scale gap. A, but a few extend out as far 
as several times the Debye energy. The only upper limit 
for CFs is the band edge, but one may anticipate that 
most states will pile up just above some characteristic 
energy scale Eco (see right-hand inset of Fig. [1]). 

The AA theory predicts the energy dependence for the 
DOS of correlated 3D metals to be for low energy, 

E, measured relative to /i^, and this dependence was con- 
firmed by several tunneling studies [31. 13. fiol. [lH . Photoe- 
mission spectroscopy (PES) showed a depressed spectral 
weight in metallic perovskites[3|] near /i^?, and recently 
the ^/E dependence was observed [lij in PES spectra of 
Sr2FeMo06. Our measurements on bilayer manganite 
single crystals, La2-2xSri+2xMn207 (LSMO) for x=0.36, 
complement all these studies by also addressing conserva- 
tion of states, a fundamental feature of electron-electron 
correlations, which has not been addressed in the litera- 
ture to date. 

Certain properties of LSMO make it ideal to study 
CPs. Low-temperature LSMO conductivity data (see 
Fig. [1]) exhibit a a/T dependence that is consistent |13j 
with AA theory. The bilayered LSMO is a strongly 
correlated system, which is readily cleaved along the 
crystal a6-plane. This enabled detailed Fermi surface 
studies using angle-resolved photoemission (ARPES), 
which show 114 Il5j renormalized Fermi velocities, vp, 



FIG. 1: (Color online) Bulk a6-plane conductivity and magne- 
tization (in a field of 1000 Oe) versus temperature are shown, 
arrows assign the scales to the curves. Right inset: Cartoon 
illustrating the expected Colomb interactions effect in 3D tun- 
neling DOS. Equal areas with vertical and horizontal lines the 
state conservation, graph centered at Fermi energy. Left inset: 
The junction geometry. 



of ~1— 3x10^ cm/s and electron scattering times, r, of 
~5— 7 fs. Such small values define a poor metal, and one 
might expect them to be due to CPs in the DOS, v^iiy)- 

The effects of electron correlations on z^(E), are 
directly observable in the tunneling conductance, 
G(V)~i^(E = eV)(xV^. A systematic analysis of the 
tunneling data yielded the values of Eco, and an observa- 
tion of a pileup of states within 2—3 times a characteristic 
energy Eco in analogy with superconductivity. Extracted 
values of vp and r agree with the published ARPES re- 
sults, and they are consistent with an intuitive model 
that estimates Eco~^/t, effectively verifying that model. 

Our tunneling data were taken on cleaved crystals of 
Lai.28Sri.72Mn207 (i.e., x=0.36) at low temperatures, 
where they are metallic. Crystals were melt- grown [l6| 
in an optical image furnace. The LSMO (x=0.36) crys- 
tals were characterised by measuring the temperature- 
dependent magnetization, M(T), as displayed in Fig. [T] 
(conductivity data are from a different crystal). Crystals 
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FIG. 2: (Color online) (a) Set of G(V) curves at 4.2 K, top- 
most dataset overlayed with its Go(V) fit. Zeroes for the 
curves indicated by bars on the right, in the same sequence. 
Zero-bias resistances of all junctions from bottom up were, in 
MQ: 80, 23, 2.6, 2.6, 1.8, 0.78, 0.46. The middle dataset is un- 
sealed (values in /iS). (b) Tunneling G(V) curves plotted ver- 
sus V^. Bottom 3 curves: The 4.2 K data to show nearly lin- 
ear dependence. Upper 2 curves: A 58 K curve scaled on top 
of a 4.2 K curve (large datapoints). The feature is sub- 
stantially thermally smeared in the 58 K curve, as predicted 
by the theory (solid black lines for 4.2 and 58 K calculated 
from Eqs. (29) and (30) of Ref. pj]). The good agreement of 
model and data reinforces our conclusion that the zero-bias 
anomaly originates from electron-electron interactions . 



(~lx 1x0.5 mm^) were cleaved in air along the easy ab- 
plane, immediately mounted in the cryostat and cooled 
down to 4.2 K. Then the gold tip was brought into con- 
tact with the cleaved c-axis-normal crystal face to cre- 
ate junctions. Details of the apparatus are described 
elsewhere [17]. 

For coherent tunneling, the reduced density of states, 
z^(E)/i/o is ~G(V)/Go(V) in the low temperature limit, 
where Go(V) is the barrier conductance that is indepen- 
dent of interaction effects and is approximately constant 
for barrier height, ^>>V. In Fig. [2^ we present low-bias, 
G(V)=dI/dV, for several junctions measured at 4.2 K. 



The data are scaled and offset (bars on the right show ze- 
ros for each curve) to show the excellent reproducibility of 
the zero-bias cusp feature across two orders of magnitude 
of junction resistances. These curves have the approxi- 
mate shape expected from the DOS effect of electron- 
electron interactions, iy/uo{E) ex a/E, as can be seen in 
the plot of Fig. [2)3 versus that shows an extended 
linear region from ~2— 10 mV. However, our G(V) data 
represent the DOS convoluted with the barrier conduc- 
tance Go(V), which must be determined. A parabolic 
approximation to Go(V) (see procedure below) is shown 
for the topmost G(V) dataset of Fig. [2^. 

Comparing our I(V) data with the WKB approxima- 
tion allows us to prove that tunneling is their origin, 
to determine the barrier height and width, and to jus- 
tify the parabolic approximation to Go(V). Freeland et 
al.[20| discovered that a non- ferromagnetic, insulating bi- 
layer occurs naturally on the surface of LSMO (x=0.40) 
crystals and the extended I(V)'s shown in Fig. [3] for 
x=0.36 closely follow the earlier data. The reasonable 
agreement over four orders-of-magnitude of current and 
three orders-of-magnitude of nominal junction resistance 
shown in Fig. [3] confirms that tunneling is the primary 
conductance mechanism. Using the collected results of 
over 50 junctions we determine that <l>=280— 380 meV 
and to =1.1 — 1.8 nm[2l[. The values of to reasonably 
approximate the full first surface bilayer thickness. Im- 
portantly, the excellent agreement with WKB tunneling 
model implies a high quality, uniform tunnel barrier. 

However, the data over the full range of V necessar- 
ily include the higher-voltage effects of barrier round- 
ing by image forces, barrier changes due to electrostric- 
tion, etc. that cannot be captured by our simple two- 
parameter WKB calculation. Thus we approximate the 
WKB model at V<<l> with a parabolic correction to 
a constant conductance, that is a current- volt age form 
Ifit(V) = a(V-Vo) + b(V-Vo)^ + aVo+bVg, where Vo 
is the voltage shift of the parabola to account for a slight 
barrier asymmetry [22|. This approximation is shown by 
the dashed line in Fig. [3] and is within 3% precision for 
V<0.7<l>/e, imposing a justifiable upper limit of ~230 mV 
for use of Ifit [23[ . 

We have fit the Ifit to our data over a variety of volt- 
age intervals within the range 0—230 mV and found that 
the Ifit is essentially independent of the interval choice, 
when the limits lie within 130—230 mV range. This is 
demonstrated in Fig. [4^, where we highlighted the differ- 
ence between I(V) data and its Ifit's using the deviation 
curves, Idata/Ifit — 1- The spread of the deviation curves 
is small and not systematic with the interval. 

In the top part of Fig. [4)3 we present a series of 
G(V)/Go(V)=i^(E)/z/o curves, and at the bottom of 
Fig. [Hd we show the integral of the {v — i'o)/^o corre- 
sponding to the top curves. Depletion of states occurs in 
the rising part and pileup in the falling part of the inte- 
gral curves. The curves show that states depleted near 
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FIG. 3: (Color online) Demonstration of tunneling: Six un- 
sealed tunneling I(V) curves (datapoints) on a log- log scale, 
with square-barrier WKB fits (solid black lines) of coher- 
ent tunneling across a square barrier. Resulting fit param- 
eters paired as to [nm], $ [me V]: #1: 1.55, 302; #2: 1.53; 325; 
#3: 1.73, 325; #4: 1.64, 330; #5: 1.67; 308; #6: 1.73, 325. 
Dashed line is a low- voltage Ifit curve fitting the WKB model 
#3. Inset: Linear plot of the top- most I(V) curve from the 
main plot in the same units. 



jiiF are recovered by ~200 meV. However, we should con- 
sider only data below 130 meV (lower limit of fit inter- 
vals) as a reasonably reliable representation of the state 
recovery. The pileup of states in u/vq plots between Eco 
and 130 meV accounts for ~80±20% of the states de- 
pleted below Eco, with the remainder presumably extend- 
ing out to the band edge as anticipated from the theory. 
The high-voltage noise, inherent to tunneling, could lead 
to the apparent asymmetry in the curves, but we have 
not detected a systematic asymmetry in the Eco values. 
From data shown in Fig. [4)3 we determine experimen- 
tally the characteristic energy scale, Eco, for the CI's in 
LSMO (x=0.36) and demonstrate that a most states pile 
up within a few Eco, in close analogy to the well-known 
superconducting case. 

The theories of interactions (AA) and quantum 
interference [isj in quasi 2D metals in the low temper- 
ature limit, allow us to estimate the electron scatter- 
ing time, r, in LSMO. The r can be found using Eqs. 
(29-31) of Ref. [T^ as 
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where a is the tight-binding coefficient for c-axis trans- 
port, and £f refers to a6-plane Fermi energy. The last 
term in Eq. ([T]) is the square of the slopes of the linear 
regions of the plot v/vq vs. a/E, see Fig. St. The pa- 
rameters Sft/H and ar/h were determined by quantum 




IOOO0— 



500 



3 
> 



-200 -100 100 200 

V(mV) 




{b)Top 




200 



-100 100 

E=eV (meV) 



200 



> 1 



ro.5 





1 1 1 1 1 1 1 1 1 1 1 1 
slopes 
(v-Vo)/vo 




-(c) 


1 , , , /tIx , , , I , 





-10 







WlE=eVI (VmeV) 



10 



FIG. 4: (Color online) Normalization procedure and results: 
(a) Example of low- voltage I(V) data (circles and line) with 
Ifit curves (10 superimposed lines), and the corresponding 
I(V)/Ifit(V)-l deviation curves. Arrows assign curves to their 
scales. Fit intervals use a fixed upper limit 230 mV, and a 
lower limit varying from 100 to 200 mV. (b) Top (left axis): 
Set of normalized conductance curves, z^/i^o, offset by inte- 
ger multiples of 0.25 for clarity. The cross-over energies, 
Eco = eVco, are marked by circles. Bottom (right axis): In- 
tegrated (y — zyo)/z^o, in arbitrary units, to illustrate how the 
pileup of states compensates the depletion near Fermi level 
(dashed Une denotes zero), (c) A single z//z/o curve plotted 
vs. VE = eV. Slopes of the linear sections highlighted by 
straight lines. 



interference [2^ as ~2.9 and ~0.15, respectively, for the 
bilayered LSMO (x=0.40). The same experiment gave an 
estimate of the mean free path £:^1.4 nm. We have an- 
alyzed ~20 junctions using the normalization procedure 
described above and Eq. ([T]). The resulting values of r 
exhibited a Gaussian distribution, with r=16±3 fs. 

Using this result we can test another, intuitive ap- 
proach to CI's. If one considers the time limits on the 
diffusive motion of electrons in the same quasi 2D model 
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as used above 25^ the result for small E is 



(2) 



where V= ^/T^'^aiPc quasi 2D electron diffusion co- 

efficient. The crossover points Eco (see Fig. [1] inset and 
Fig. [4)3) between depleted states near (ip and the pile- up 
at higher energies is defined by v=vq- Then using r from 
Eq. ([1]) we determine Eco=^/t from Eq. ([2]) and com- 
pare it with the experimental values shown in Fig [4)3. 
The AA theory above predicts Eco=41±8 meV, which 
matches with the Eco values obtained from Fig. [4)3, i.e., 
33—55 meV. This verifies the intuitive formula ([2]). 

The theory of interactions has been developed as a 
perturbation theory for small corrections to the DOS, 
in the limit of ^/Af>>1, where ^=vft is the mean free 
path, and Ap is Fermi wavelength. The DOS effect we 
present in Fig. [H is not so small, and the stated assump- 
tion is barely satisfied, as ^/ \y=Syt /h^2.^. The val- 
ues of T determined by our analysis, therefore, are not 
rigorously correct. Nonetheless, they are in a reason- 
able agreement with recent ARPES results [13, EH. Us- 
ing r=16 fs and the Ic^lA nm, we find a Fermi velocity, 
vp = ^/r ~ 1x10^ cm/s. This is close to the published 
ARPES results, vf=1-3x10^ cm/s. Note that ARPES 
experiments measure the Fermi velocity along selected 
directions in the /c-space, whereas our tunneling experi- 
ment measures vp averaged over the /c-states, weighted 
by their contribution to tunneling. 

In summary, we find that ~80% of states near (ip de- 
pleted by Coulomb interactions are found in the range 
of ^40—130 meV from fip. This presents an analogy 
with superconductivity, where a large fraction of depleted 
states reappear within two or three times the character- 
istic energy. However, the characteristic energies, Eco, 
are not universal, e.g., they depend on r in LSMO and 
on the energy gap. A, in superconductors. Our anal- 
ysis yields elastic scattering time, r = 16 ± 3 fs, and 
Fermi velocity, vp = ^/r ~ 1x10^ cm/s, that are not too 
far from recently published ARPES results. We accom- 
plished this by a consistent procedure to normalize the 
tunneling conductance. This appears to be the first ex- 
periment to address an essential implication of the AA 
theory of CPs, the conservation of states. 
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